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Abstract 

The QCD corrections for the box diagrams are revisited for the 
case of a heavy top quark with mt = lYAGeV. We resolve first a long- 
standing discrepancy between two methods of calculation by showing 
that they give the same results when the threshold factors are treated 
correctly. Using this observation we refine our earlier results and de- 
rive formulae valid for the K- and S-meson systems. Our formulae 
are given in terms of integrals to be evaluated numerically, as well 
as approximate analytical formulae. These calculations include the 
evolution above mw which has been neglected by other authors. 
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1 Introduction 



Perturbative quantum chromo dynamic (QCD) corrections to electroweak 
processes, dominated by short distance physics, were calculated in the leading 
logarithmic approximation (LLA) long time ago The pioneering calcu- 
lations studied the enhancement of the AI = 1/2-amplitude in AS* = 1 
processes and the origin of direct CP violation. 

Subsequently, QCD corrections in the LLA were computed for AS = 2 
processes [0-0] assuming that all the quarks in the loops are much lighter 
than the W bosons (m^ <C mw for all q). In course of time, the gradual 
strengthening of the lower limit on the mass rrit of the t quark lead to the 
realization that rrit may indeed be comparable to mw or even larger. This 
motivated several papers which generalize the calculations for a heavy top 
quark 0-0- Some of these articles 0-0 simpified the calculation by ne- 
glecting the evolution of the Wilson coefficients in the interval mw mt- 
This approximation was not made in However, since mt was unknown, 
the numerical significance of this approximation could not be assessed prop- 
erly. 

In the meanwhile, the production of the top quark has been reported by 
two Fermilab groups. The CDF collaboration [|T^] finds a top quark with 
a mass mt = 174 ± 10 ± ^^2^61^ and the DO collaboration |ll[ with a mass 



mt = 180 ± 12GeV . These values are in good agreement with the precision 
measurements at LEP-I where it was established |13| (albeit indirectly) 
that mt = 169 ± ± ^J^GeV. All this forces us to accept that the mass 
difference between the top quark and the W-boson is indeed significant. It 
is therefore worthwhile to revisit the calculations of and compare them 
with other results in order to estimate the numerical significance of the 

evolution beyond mw- We find that this evolution is numerically significant 
for mt ~ llAGeV although it does not change the result drastically. 

The first isssue concerns the method of calculation to be adopted in this 
article. There are two methods : the first one is diagramatic where one loop 
QCD corrections are computed on top of the weak box diagram. The lowest 
order corrections thus obtained, are summed up using the renormalisation 
group (RG). The resulting Wilson coefficients are dependent on the momen- 
tum of the electroweak loop (p^) which are then integrated over 0, ^, |]. 
The second method first introduced in H] computes the QCD corrections in 
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a series of effective field theories obtained by successively integrating out the 
heavy degrees of freedom 0, ||, ^ |^ . The two methods have a long-standing 
and unsettled issue which was raised in ref. 0. To be specific, even in case 
of rrit <C myy, the results of (which heavily uses the rather involved tech- 
niques of operator mixing) and ref. did not agree. Indeed it was noted in 
ref. that although two of the three QCD correction factors(?7i, 172) com- 
puted by them agree with ref. 0] after appropriate simplifications, the third 
(773) did not (see footnote 14 of ref. 0). To the best of our knowledge this 
discrepancy has not been clarified in the subsequent literature. Since the 
disagreement was precisely in a term where operator mixings played an im- 
portant role in ref. 0, a critical reader may also question the procedure of 
ref. [Q, which did not explicitly use operator mixing. This critique would 
also apply to ref. [|5| where the techniques of Novikov, Shifman, Vainshtein, 
Zakharov ^ and Visotskii [^] (NSVZV) was generalized for nit ^ ttlw. 

In this paper we explicitly demonstrate that the apparent disagreement 
between ref. and 0] is a consequence of certain simplifying assumptions in 
ref. More specifically in ref. [^] certain Wilson coefficients were evolved 
over large momentum intervals which involve several thresholds. Some re- 
marks were made regarding the ambiguities in determining the number of 
active quark flavours (n/) in these intervals, but for the purpose of actual 
computation, Uf was kept fixed. We find that by changing n/ at each new 
quark threshold and taking the matching conditions at each threshold |Q into 
account a more rigorous formula can be derived from the procedure of ref. 
[^] . This leads to an agreement — not only numerically but also algebraically 
— with ref. @. 

Although the case rrit <C mw is now of academic interest only , we con- 
sider the above demonstration important because it illustrates that both 
methods of calculation (GW and NSVZV and its generalization 0) are equiv- 
alent. The apparent disagreement comes from the approximate treatment of 
the threshold factors as discussed above. In the present article, we remove 
these approximations and sharpen the formulae of ref. [^. In section 3 we 
present explicit and improved formulae which are in one-to-one correspon- 
dence to our previous work The formulae by themselves answer several 
questions raised about our earlier work (see in particular |^). 

The main new result of ref. was to demonstrate that a particular 
Wilson coefficient computed in 0, |] does not contribute if >> mw The 
remaining Wilson coefficients were calculated. The dependent Wilson 
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coefficients as discussed above were integrated numerically over the entire 
permissible range of p^. In this paper we incorporate the threshold factors 
discussed earlier, into these integrations and obtain very accurate results. 

However, from the final results thus obtained the numerical significance 
of the evolution of as{p^) over a specific interval {my/ — > rrit) cannot be 
isolated automatically. This can be done by separating the regions of the 
numerical integration and setting a^d?^) = ^^(m^) for > m^. Alterna- 
tively, in order to separate the evolution of this specific interval, we introduce 
a reasonable approximation, which allows us to perform the p^-integration 
analytically. This way we isolate the interval {mw — > m^), where we can 
choose as(p^) as either running or as a constant with as(p^) = as(m^) for 
> royy . We thus obtain for the ffist time a set of analytical formulae which 
explicitly exhibits the effects of the evolution above rriw and readily reduces 
to the results of [7|, ^ once these evolutions are neglected. However, as 
we summarize below in various tables, the evolution above mvi/, though not 
completely negligible, does not change the results drastically. This is partly 
due to the fact that various corrections cancel each other. Finally, the for- 
malism of this article is very transparent and general so that it can be easily 
extended to four generations, where higher scales become relevant. In such 
cases the evolution above mw can in principle be crucially important. It 
has been shown that models with naturally heavy Majorana [II^ or Dirac 



15] neutrinos belonging to the fourth generation can be constructed. Such 



models are ,therefore, neither unnatural nor in conffict with the neutrino 



counting at LEP|2|, |T3 



In ref. p next to leading order calculations for one of the QCD factors, 
r]2 in the notation of P], have been performed. These attempts are important 
since they have the potential of settling important theoretical issues, like the 
independence of the final result from the scale at which the initial conditions 
for the Wilson coefficients have been imposed or from the definitions of the 
heavy top masses, etc. Such refinements however can attain their full poten- 
tial only if the corrections in the LLA are known as accurately as possible. 
In this paper we find that in some cases the evolution above mw is numeri- 
cally of the same order as the nonleading corrections and reqires, therefore, a 
closer scrutiny (see below). This explains, in part, the fact that the authors 
in 1^ obtained for 772 the same result. For the case > mw there is for 773 
only the calculation in ||^ and the present one. 

The paper is organized as follows. In section 2 we briefly review the 
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well known basic formulae for AF = 2 processes relevant for — and 
— mixings in order to set up the notation and make the article self- 
contained. We then review the computation of the QCD corrections in the 
hypothetical case <C following the NSVZV method and demonstrate 
the equivalence with GW. This is reached by treating the threshold factors 
and the matching conditions carefully. Using the experience gained in sec- 
tion 2, we refine the results of ref. and present them as specific integrals. 
This is done explicitly for each integral of the effective weak Hamiltonian 
and summarized in section 3. In addition, section 3 includes analytic for- 
mulae obtained after integrating the terms, in a certain approximation, in 
order to separate the evolution above mw and compare it with other works 
[@] ^ [§!• The formulae look long and complicated, but are obtained by a 
straight-forward treatment of thresholds. In section 4 we evaluate the results 
numerically arriving at final values for the QCD factors and compare them 
with those in other articles. 



2 QCD Corrections in the NSVZV approach 

In this section we consider the hypothetical case my/ ^ mq (for all q) in 
order to connect with the early work on this subject and in particular to 
demonstrate the equivalence between the works of Oilman and Wise and 
Visotskii If. 

In this limit {my/ ^ mq) the effective Hamiltonian without QCD correc- 
tions is obtained from the box diagram (Fig. 1) with two W boson exchanges 
and attains the form |[T6[| : 

(72 ^2 

^ ^^(^7.(1-75)^)^ E A.A,if,^/(^,j) (1) 

■^"^ i,j=c,t 

where the GIM subtraction is explicitly introduced with mu ~ 0. We follow 
the standard notation Aj = V*iVis, where the V"'s are the elements of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix. The dimensionless integrals 
I{mf^m'^) are over Euclidean loop momentum. 
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^ ^' i (p2 + m?)(p2 + m2)(p2 + ^2^)2 

These integrals are convergent so that the dominant contributions come 
from p'^ ^ mf, as well as p"^ ~ rn^^. Due to the GIM subtraction in Eq. (|l|) 
the contribution form the region cancels out. Hence it is justified 

to replace the W propagators by m^. Now Hefj is implicitly the time- 
ordered product of two AS = 1, local four fermion operators. This of course 
can be seen most transparently in coordinate space (see Eqs. 5a-5c of ref. 
[§]). These operators can be written as linear combinations of the well known 
colour symmetric and antisymmetric operators O^. The QCD corrections for 
these multiplicatively renormalisable operators are well known Finally 
the QCD corrected AS* = 1 operators (neglecting mixings with penguins) 
are joined again by Wick's rule. In momentum space this corresponds to the 
replacement of the W boson propagators in Eq. (|^) by and multiplying 
each integral of Eq.( [I|) by a factor [0, ^ 

Qi{mw,P,ni) = -{ r^)'*"!'-! T^)^"'' + 



2 ttni(p^) ani(p^) 2 an^{p 



,2^ 



I 6(11) 



(3) 

where a+ = 2, a_ = —4 are the anomalous dimensions of in units of 
a„^/(27r), a„j is the strong coupling constant and b{ni) = 11 — |ni depends 
on the number of massless quarks in the interval mw — ^ p"^- We emphasize 
that since p"^ is variable Ui changes as p"^ crosses various quark thresholds. 

The two momentum scales in Eq. (^ are worth noting. It was emphasized 
in ref. that for Green functions involving loops, provides an important 
scale and the emergence of momentum (p^) dependent Wilson coefficients 
was explicitly demonstrated ( see in particular the third article of ref. 
Moreover, as pointed out in ref. 0, ^, and directly shown in ref. [Q, the 
0{a) terms in the expansion of Eq. (^ can be directly identified with the 
diagrams of figure 1 (see Eqs. 20-2 lb of ref. p). 

Once p^ becomes larger than the masses and momenta on the external 
legs, the effective four fermion operator {dj^i^l — 75)5)^ is generated. This 
property is already implicit in Eq.( |l|), which is derived in the approximation 
of neglecting the external momenta and it also holds in the realistic case 
rrit > my/. The evolution of the four fermion operator ((^7^(1 — 75)5)^, which 
is colour symmetric and has the same anomalous dimension as 0"*", from p^ 
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down to /i^ (a scale where perturbative QCD presumably breaks down) yields 
a factor 

Q2(p^/x^n2) = (^^)^ (4) 

"nal/^ ) 

which multiplies the integrands of Eq. (1). The 0{a) term in the expansion 
of Eq. (^) can be explicitly checked from fig. 2 as was done in ref. (see 
Eq. (12) of and the discussion following it). It is very important to note 
here that the number of active flavours (n2) in Eq. (^ (which is not fixed as 
yet) can be different from that in Eq. (3) since the momentum intervals in 
the two cases are different. 

Finally each of the heavy quarks in the internal lines are off mass-shell 
during the integration and can be dressed with gluons as in figure 3. This 
corresponds to replacing each heavy quark mass in Eq. ([^) by 0, ^ 

mlip') = mj(mj)(l + 6(n3)^^^ln4)-'/'^"^^ 



q 



a„3 (p 



,2^ 



which is the well known relation between running quark masses at scales 
and m^. The right-hand side of Eq. (^ was not shown explicitly in ref. 
[0, ^, 1^, but was evidently defined correctly. This caused some criticism in 
the subsequent literature p|, which is not justified. 

Following our above observations we have kept the number of flavours 
(n3 in Eq. (^) arbitrary and in principle different from those in Eqs. (0,§). 
It is rather obvious that this correction is also independent of the condition 
m^r ^ mf and continues to hold in the realistic case m,t > mw- 

The removal of the W boson propagators and GIM subtraction leads to 
a generic contribution to the effective Hamiltonian defined in Eq.( 1) (in the 
absence of QCD corrections) 

^^^^ '•^^ ml. J {p^ + mf){p^ + m])- ^ ^ 

When two quarks on the internal lines of the box diagram are identical 
(i = j = c or t, henceforth referred to as the cc or tt graph) the integral is 
dominated by contributions from ^ mf [El ffl. On the other hand, for the 
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diagram with one internal c quark and one t quark (the ct graph), the entire 
region < < rnf contributes as can be seen from the following form 



Heff{c,t) 



mim^ f"^l dp- 



x.ln4 (7) 



where Xq = m^/m^. 

Due to QCD corrections the above formula should be modified to 

1 /• dp^mf(p^)m|(p^)gl(m^,p^nl)g2(p^/^^n2) 
~ ml, J (p2 + m2(p2))(p2 + rn]{p^)) ^ ^ 

where the corrections are defined in Eqs. A very accurate analytical 

result can be obtained by the following approximation. It is assumed that 
dependent QCD corrections cannot drastically alter the important regions 
of momentum integration as noted after Eq. (^. For the i = j case, this 
corresponds to evaluating the corrections (Eqs. at p"^ = m]{mf) and 

one obtains from Eq. (§) 

Heff{i,i) = x^Ql(m^,m^,^^l)Q2('^^At^"'2) (9) 

For notational convenience we shall denote henceforth mj(m^) = mj. 

In ref. @], n2 in Eq. (^ was assumed to be 3 for numerical calculations. 
While this is justified for the cc graph the same cannot be said about the 
tt graph. Similarly rii was chosen to be 4 which is not rigorously justified 
either for the cc or for the tt graph. For the cc graph, for example, one can 
replace a typical term in Qi of Eq. (H) by 

^ a„^(m^) y/b(ni) ^ ^ (yejfnw) y/bje) ^ Q^s (^t ) y/b(5) ^ Qj^i'mf) -^g/bi^) /^^v 

where the changes in the number of flavours at different thresholds have 
been explicitly included. It is further assumed that the matching conditions 
as defined in 0] are also obeyed. With this modification Eq. (^ reproduces 
the results of ref. , which was already noted in ref. |0 . A similar agreement 
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for the tt graph can be obtained analogously by modifying the Q2 term in 
Eq. i). 

The situation is somewhat more complicated for the ct graph. Using Eq. 
(1^, we can write the integral of Eq. (P) in the following approximate form 

Hf^^'ict) = r' ^m^(/)QlK^,p^nOQ2b^/i^n2) (11) 

When we substitute in Eqs. (^1)-® the number of flavours ?^i,2,3 = 4, 3, 4, 
respectively, and factorize 

and 



2^ 



)(-7-|f) 



which roughly corresponds to the approximation that above rric the number 
of active quark flavour is fixed at 4, the integration in Eq. (|Tl]) can be done 
easily and 773 of ref. is rederived. The resulting formula unfortunately 
disagrees with ref. 0. This is the disagreement already mentioned and will 
be resolved below. 

In order to treat the number of active flavours consistently, we note that 
a typical term of Eq.( 0) is of the following form: 



Jmi p2 ansiml) ttnal/W^) "^("^Vf) 

Values for rii, n2, and 113 can be assigned consistently by splitting the 
integral into two pieces < < ml and ml < p^ < ml . Now each 
momentum interval can be assigned a well defined number of active flavours. 
This yields 



, a3(mg) (2/b(3)) / O^ei^t) \(a/b(6)) 

f /M';«|)xa/6(5)/M^Na/M^ ^ . Mp^x (8+a+2)/fe(4) 



'a5{ml) ce4{ml) Jml p^ ai{m'l) 



ai{m^^} OLb{mt} Jm^ p^ a^{mi 
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where a= -4, 2 or 8 and replacements similar to Eq.( 10) have been made. 
Using the general formula 



J{ml, ml, x) 



p2 



h{n) 

477 



a„(mi) In^ 



477 



Mn{ml) , 

6(ra)a„(m|)(l - x) an{mj) ' 



(14) 



the integrals in Eq.(|T3D can be done analytically. 



In simplifying the results, the matching conditions (e. g. a4(m^ 



[§]) should be used in particular for the anifri^) in the overall factor in Eq. 
(0). This leads to, e.g. 



«5(m?) = (^)("^(-^^ 



Using all of the above one finally obtains 

_ , a3(^c) N2/fe(3)/ M^t) Na/fe(6) 



)a4{ml) 



477 



(15) 



a3(/i2 



1 - 



a4(m2)(6(4) - 10 - a) 

2 
b 
2 



&(4)-10-a \ ^ a4(r»g) ^io/fc(4)_i^ a5(^?) ^-a/b(5) 
6(5) — 10 — ay 0:4(771; 



'05(777^) 



' «5(w^f) X -a/fe(5) / «4(mg) _„/fe(4) 



'0:4(7^^) 



+ 



^ &(4) - 10-a ^ a4(m|)_^io/f,(4)_i "5(^^^10/6(5) - 
6(5) — 10 — a 0:4(7^7^) 05(7715) 



(16) 



In this notation Eq. ([TT|) can be rewritten as 



Heffic, t) 



(17) 



which is algebraically identical to Eq. (23) of ref. In summary we reiterate 
that the effects of the mixing of the operators O7, Os in the coordinate space 
as studied in ref. [§, can be equivalently treated in this momentum space 
calculation by handling the integration properly (notice the mixing of the 
anomalous dimensions of various operators in the integrands of Eq. (13)). 
Consequently, one can use the methods of GW or NSVZV and obtain the 
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same results, provided that threshold effects are handled appropriately. In 
the next section we include these effects in the approach of NSVZV in order 
to refine the results of ref. for the case of rrit ^ mw 



3 QCD corrections in LLA for nit ^ t^w 

We now present the QCD correction factors for the realistic case rrit 3> mw 
which is a refined version of the calculations of ref. along the lines devel- 
oped in section |^. In this calculation diagrams with unphysical Higgs scalar 
exchanges are also important. In the following we denote by the subscripts 
WW, WH and HH the contributions of diagrams with two W boson internal 
lines, one W and one Higgs boson internal lines and two Higgs boson internal 
lines respectively. For the cc diagram the approximation mw ^ is still 
justified. Thus the results of ref. ^ are still valid. 

For the ct graph the terms in Eq. (0) can be regrouped in the following 

way 

H^w{c,t) = [l{0,0)-I{0,ml)] (18) 
H'^^ic t) = [I{ml ml) - /(O, m?)] (19) 

In Eq. (|18D which is independent of mj, the W boson can be considered again 
as heavy. However, since there is only one GIM subtraction, both W boson 
propagators cannot be removed directly. Instead the heavy W limit can be 
realized by following the steps leading to Eq. (0), i. e. 



dp^ 1 f dp 



ijp' + m'DijP' + m\rY mw — m^ J p"^ + p"^ + m^ p"^ + 



- dp^ 



f^W "^W "''"c "^w 



where he first intgral recieves contributions from < < Tn^r and the 
second one from p^ ^ m\r'^ ^ Thus we obtain the QCD corrected form 



WW 



dp^ 
p2 



-m; 



'{p')Qi{p')Q2{p') - ml{m'w)Q2{m'w) (20) 
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This was already derived in ref. where the number of flavours, however, 
was treated in an approximate way as discussed in the last section. Following 
section ^(see Eqs. 11 - 17) a rigorous result can now be obtained from Eq. 



(pOD-The integral is the same as Eq. ([TT|) with replaced by mw and 
some obvious readjustments of the number of flavours. After these minor 
modifications the result agrees completely with the corresponding formula of 
ref. IP] which , however, neglects the second term in Eq.( 20) (see below for 
further discussions on this point). 

For Eq. (|19]), however, the heavy W limit is no longer justified. In fact 
as noted in ref. P, the corrections in Eq. (^) do not arise any more, because 
the diagrams with the gluon connecting an external with an internal quark 
line decouple in this case. Then Eq. (|T^ with QCD corrections becomes 



Here only the leading order term in is retained. For accurate result the 
integration in Eq.( 21) is done numerically considering the region p"^ > fi^ 
only. In course of the integration the number of active flavours in the QCD 
corrections should be adjusted according to the value of as discussed in 
section 2. 

An analytical formula can be easily obtained by the replacement mf{p'^) ~ 
mf in Eq.( 21). This corresponds to the approximation {mf{p'^) —Tnf)/{p'^ + 
m^) -C 1. Since for the above integral the dominant contribution comes from 
iTTw ^ ^ ^1-1 ^his approximation is reasonable. Now splitting the inte- 
grand such that different terms receive contributions from different regions 
of p2 I, I, H we obtain 



H\y\y{,c, t) 



dp^ml^m^ 1 1 



rni - {p^ + mf) {p^ + tu^) {p^ + m^Y 

■j2/6(n2)^ '-^"■3 (P^) N8/fe(n3) 

1 1 n 

[— — J(m?,m^,-)-l] (22) 

Xt — I O5 



As already discussed (see e. g. Eq. (0)) the second term in the square bracket 
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of (^) can be obtained by evaluating the QCD corrections at p ~ Tn^r■ 
The first term, however, should be evaluated as in Eqs.7, 13 and 14. This 
procedure leads to the approximate result of eq. (P^). 

In ref. QCD corrections in the interval were neglected. In 

this approximation (a(p^) ~ a{m'^) for > m^), J in Eq. ( p2D reduces to 



In Xt and the result of ref. is regained. The evolution above the scale mw 
is given by J. However, J multiplies a term suppressed by Xt — 1 and does 
not affect the numerical result appreciably for rrit = VJAGeV . It should be 
noted that the new correction is not an overall multiplicative factor. This 
feature can be seen in all the subsequent formulae. 

In ref. it was suggested that QCD corrections to the ct graph can be 
computed by substituting mw ~ rrit in the corresponding formula of ref. [Q]. 
In view of Eq.( 20) and the discussions following it, it is clear that this is 
justified for the part of H^ff given in Eq. (|TBp, apart from the neglect of the 
second term in Eq.( 20). We have checked that this corresponds to an error 
of 6 %.This prescription is obviously not justified for the piece in Eq.(19). 
However, the contribution of Eq. (19) is numerically much smaller than that 
of Eq. (18). A posteriori the approximation of appears to be reasonable. 

Following the same procedure as above one obtains (keeping only the 
leading order term in ml) for HH and WH exchanges in the ct diagram: 



p2 



(p2 + mlij)^))^^^ + rriy^Y 



(23) 



and 



1 



(p2 + (p2))(p2 _^ ^2^j 



(24) 



The approximation for performing the integrals analytically has already 
been described . Introducing this we arrive at the results 

Am^ J ansimi) oin^A^t) "^^(/i^) 
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+ 



XqXt / Q;3(m^) x 2/6(3)/ 0^4(^6) -.10/6(4) / Q^5("^h/) Nio/b(5) 



1 .a5(m^) 



- 1 ar,{m1) 



\8/b(5) 



and 



1 



1 



{ml - m'^){p^ + mi2)(p2 + ^2^) 
23:cXf . 0:3(772^) 2/b(3) , Q;4(mj,) xio/b(4) . ^sl'^^ty) Nio/b(5) 

r/ Q5("^vt/) x8/fe(5) _ 1 / Q;5(m^) g/j,(5) 2 2 J^- 



(26) 



As has already been stated, they give good estimates, but for very accurate 
computation it is desirable to integrate the exact equations numerically. The 
QCD corrections from my/ are contained in the square brackets. The 

three overall factors give the QCD evolution from /i fnyy, and agrees with 
ref. 

Finally, we consider the tt graphs. For reasons discussed earlier, the 
correction in Eq. (^) are no longer relevant. Including the other two types 
of corrections we obtain for the WW, HH and WH exchanges 

(27) 



Hwwit.t) = ml/ j dp^m^{p^){ 



(p2 + m1{p'^)Y{p'^ + my/ 
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,,.4 . ^, (^2 

HnHit, t) 



/'^p2/ '^n2(P ) x2/b(n2)/ "^"S (P ) \16/&(n3) 



(28) 



\2/b(no.)/ ^nz 



vl6/6(n3) 



(29) 



(p2 + (p^))2(p2 + m^y) 



The two types of QCD corrections arise from the running top quark mass and 
the diagrams in fig. 2. They produce the terms containing 6(713) and 6(^2), 
respectively. Again the numbers of flavors should be assigned consistently at 
each momentum interval. Thus these formulae with the appropriate choice 
of the number of generations n2,n3 and fi can be used for computing the 
contributions of the box diagrams for K- and B-mesons. 

Introducing again the approximation described after Eq. (21) we arrive 
at the formulae : 



[ I + I 2 ^ ], 

{p'^ + mjy {p'^ + niwy (p^ + mfj^) (p2 + ) 

_ / Q^3(^c) x2/b(3)/ ^4(m^) ^ 2/6(4) / «5(^Vk) ^ 2/6(5) 
[■ , Ctbi^D \2/6(5) _|_ 

-2 (^.^'i\^3 J(m^m^,^)], (30) 

Special care should be taken in deriving cq. (30) from eq. (27). The m^(p^) 
in Eq. (27) originates from the top propagator due to GIM subtraction. 



15 



Consequently, when we neglect the running of the top quark mass in the 
propagator, we must consistently neglect the running of the m^(p^)-terms. 
Keeping this in mind we also obtain 



{p^ + mlY {p^ + mwY {p'^ + m^y){p'^ + rrii) 

' «3(^c) 2/6(3) / «4(mg) ^2/M4) ^A^^l^ ^m^) 

" asijj'^) a4{ml) a^imf) 



i , (^5{mf) 2/6(5) , ,«5("^^ 



)2/6(5) _^ 



WJ\ie/b{5) 



4(xj - 1)2 a5(m^) 4(xt - 1)^ a5(m|) 



,,3^^^T^)"/^^'^^(-^-'v.,^)], (31) 
2{xt - ly a^imf) h 



H^nit.t) = 2mtf . , , (^^)2/Kn.)(^n3iP^)ia/Kn3) 



{p^ + {p^ + rawY 

(p2 + m2^)(p2 _^m?) 

/ M^x2/M3) .M^x 2/6(4) f^5irr>w)^2/b(5) 

^a3(/^^)^ ^«4K2)^ ^«5K)^ (xt-iy 

^(xi-l)^«5K) ^ ^ *' ^'6(5)^ 

_/«50;«w)\-2/fe(5) _ /^^sO^^^^wOn 16/6(5)1 fon\ 

^5K)^ ^sK)^ J- ^''^ 



We note that the three overall factors in Eqs. (30) - (32) give the QCD 
evolution from mw — > yU as was already noted in ref. |p, §]. The results of 
ref. P] are obtained by setting a{m^) = a{m^) and J = In xt in Eq. (30) - 
(32). 
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For B — B mixing only tt diagrams are relevant. QCD corrections can be 
obtained by carrying out the integrations in Eqs. (27)- (29) numerically. 
However, in this case ji — 0(1713) ■ Hence only > ml should be considered. 
The analytical formulae in this case are 



TT (/ /^ ~ , a5(m^y) 2/b 5 r/ Q5(^i) N2/b 5 I ^ 

a^{lJ?) [xt - 1)2 a^im^) 

(33) 



TT _ a5(m,y) 2/b(5) \^ t (^^i^t) \2/b{5) I / (^^{^w) \lQ/h{5) 

nHH[t,t) w ( 7~2^) T\2 ^^*( — T^^^ + ( — 



2Xt ,a5("^^)M6/fe(5) T/_2 _2 1 



)WM^^J(m^m^^,f)], (34) 



HwH{t, t) 



, 0:5(772^) 2/b(5) 2a:^ 



c^bii^'^) {xt - ly 

r(l+^.«5(rn|^.i6/b(5) j(2 2 18 X 
[Xt - 1) a^{mi) 6(5) 

^ Q;5(m^) ^_2/j,(5) _ ^a5(?7i|^^ 16/6(5) I ^gg^ 



0:5 (mf) Q;5("^t) 
We are now prepared to discuss the numerical results. 



4 Numerical Results and Conclusions 

The main results of this article arc QCD corrected effective Hamiltonians in 
the LLA for — K^ and B^ — B'^ systems given by the formulae (23)-(35) in 
the case rrit > mw- We present for the first time analytical formulae which 
show explicitly the effect of QCD evolutions in the interval mw — fnti 
neglected by other authors. Further we have shown the agreement between 
different calculations which apparently use different methods. 

All parameters in the formulae are defined except for the infrared scale 
/X, which is a substraction point. For physical states the low energy scale 
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must be set by the masses of the quarks and the dynamics of confinement. 
A reasonable expectation is = O^ttik) for K-mesons and /i = O(mB) for 
B-mesons. This, of course, means that the renormahzation for B-mesons 
is stopped at Mb, while for K-mesons is continued down to and below the 
rric threshold. For this reason formulae for the two cases are different: for 
B-mesons we do not cross the nic threshold. 

The /i-dependence brings in a factor (a^jl/^^))"^^^*^"^^ which should be 
cancelled, in principle, by a corresponding factor in the hadronic matrix ele- 
ment. This motivated some authors [|, ^ to factor out such a term and define 
the QCD factors without it. This does not solve the problem but pushes it to 
another part of the theory; up to now there is no explicit calculation where 
the /i-dependence of the matrix element cancels the /^-dependence from QCD. 
For the K-meson we will take Q;s(m|-) ^ 0.29 and then take the 6/25 root of 
it, which brings the /x-term closer to 1. For B-mesons we will take /i = 0(m£). 

We use the formulae of section 3 to calculate the integrals for vari- 
ous values of a heavy top quark. We present our results in several ta- 
bles. Table 1 shows the corrections for the tt-graphs of K'^ — mix- 
ing for mt = nAGeV^rric = 1.3GeV,mb = 4:A9GeV,mw = SOGeV and 
fi = 0.50GeV. Throughout the paper the Aparameter of QCD with three 
active flavours will be taken as input with a value A3 = 0.3 which is con- 
sistent with the value of A5 in [|12|. This parameter for other numbers of 
active flavours will be determined fron the input value via the matching con- 
ditions 0. We show three terms of the Hamiltonian with WW, WH and 
HH exchanges, respectively. The first column is the box diagram with only 
weak terms. The second column gives the values when QCD corrections are 
included and the integrals in Eqs (27) - (29) are evaluated numerically. As 
we discussed in the previous section, the evolution of tts(p^) for > m^r is 
now included. We notice that the change from the first to the second column 
is large, giving an 

r/2 = 0.56. 

This means that, as already noted by previous works, the strong contributions 
are sizeable and must be included. 

In the previous sections, we also demonstrated that we can calculate 
the integrals analytically by introducing an approximation. We calculated 
these terms and found out that they are in reasonable agreement with the 
numbers obtained by numerical integration. The analytical formulae are 
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particularly useful for studying the evolution above mw because the effects of 
this evolution is explicitly separated. The separation shows that the evolution 
above myy is not quite negligible when corrections to individual terms in Eqs. 
(30) - (32) are considered. For the purpose of illustration we show below the 
numerical values of the QCD corrections in Eq. (31). 



HHHit,t) ^ 0.58 



4(xt - ly 



2xt 



[0.97x4 + 1.29-1.12 —Inxt], (36) 

In Eq. (36) the pure electroweak term is kept unmodified, i.e. setting all 
numerical factors equal to one we recover the electroweak term. The factor 
0.58 is due to the QCD correction in the interval /i^ — > B "^^^ 
effects of the evolution in the interval mw — ^ are shown explicitly as a 
modification of the coefficients in the square bracket. The bulk of these effects 
comes from the running of the top mass which in this case has significantly 
larger anomalous dimension compared to the corrections in Eq. (4). The top 
quark mass correction is absent in the first term in Eq. (36) since this term 
is dominated by p"^ ~ mf (see the first step of Eq. (31)). The second term 
which is dominated by p"^ ~ (see the first step of Eq. (31)) receives the 
maximum correction (1.29) , as expected. But the overall numerical impact 
of this term is rather modest since it multiplies a small electroweak term. 
The complete results are shown in table 2. The first column gives the results 
from the full analytical formulae ( Eqs. (30) - (32)) while the second column 
is obtained by ignoring the evolutions above mw- We also performed the 
numerical integration in the latter approximation (i.e, by setting a(p^) = 
a(m^) for < p^). The results are given in the third column of table 
1 and are in fair agreement with the analytical results. The corrections for 
the Hww and Hhh terms differ from those on the second column, but the 
deviations are in opposite directions so that they partly cancel in the sum. 

In Table 3 we show the results for the ct-diagrams for mt = 174GeV. We 
note that the largest contribution comes from the H^^■^Y term (Eq. (20)) whose 
integral is dominated by the region < < m"^ . All other terms are 
smaller. We obtain the correction rj^, in agreement with previous results. 

Table 4 shows the corrections for the tt-graphs of the K-mesons for mt = 
150 and 200GeV/c^. We show again three terms of the Hamiltonian with 
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WW, WH and HH exchanges. We emphasize that the values for 772 are the 
same demonstrating, a posteriori, that a precise definiton of the running 
rrit mass is not crucial as long as it is with in the range permitted by the 



experiments 10, 11 



We did not calculate the ct-graphs for other values of rrit because the 
dominant term by far is H^^r which is independent of rrit. In summary 
we see that the QCD corrections to box diagrams of the K-system are well 
understood and give precise values. The values for the — system are 

rj2 = 0.56 and 773 = 0.36 . 

We computed the effective Hamiltonians for the B-system also for three 
masses of the top quark rrit = 150, 174 and 200GeV and for the parameters 
me = l.SGeV, mb = AA9GeV, mw = SOGeV and /i = AA9GeV. In Ta- 
ble 5 we show the three terms of the Hamiltonian for WW, WH and HH 
exchanges, separately. We note that ri2 is now larger than the K system. 

ri2 = 0.81 

and stable on various mass of m^. A large part of the change comes from the 
change of the scale fi.ln the last column of Table 5 we give the results for 
a hypothetical T quark belonging to the fourth generation. In this case 772 
changes perceptibly. Ignoring the evolution in the interval mw — > rriT one 
obtains 772 = 0.84. We have also verified that the analytical formulae Eqs. 
(33) - (35) give numbers in agreement with Table 5. 
We can compare these results with other calculations. 

1. For the K-system. The value of 7/2 is practically the same as in our 
earlier paper p. For 772 there is a next-to-leading order calculation 
which gives again the same value. We have shown in this article that 
7/2 is stable under variations of the top quark mass between 150 and 
200 GeV . For 773 the new value is in agreement with our earlier result 



2. For the B-system the dominant contribution comes from the tt-diagrams. 
The 7/2 is larger than in the i^^-mesons. This comes from the change 
in the /i-scale and the higher threshold. In comparing with ref. |^ one 
must be careful to correct for the (0772 (z^^))^''^*'^^'' which is included in 
our formulae. In ref. this factor is extracted and consequently the 
value for 7/2B is practically equal to r]2K- 
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In conclusion, the QCD factors which enter calculations with box dia- 
grams have stabilized over the past few years. Improvements through the 
inclusion of thresholds and variations on the mass of the top quarks dis- 
cussed in these papers change the values very little or not at all. The QCD 
factors are much more stable and better understood than the other param- 
eters, in particular the Cabibbo-Kobayashi-Maskawa matrix elements and 
reduced matrix elements, which enter these calculations. 
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without 


with QCD 
numerically 


with QCD 
negl. above mw 


Hww 


0.61 


0.37 


0.36 


HwH 


1.25 


0.72 


0.73 


Hhh 


0.72 


0.36 


0.42 


sum 


2.58 


1.45 


1.51 


r]2 = 




0.56 


0.59 



Table 1: QCD corrections for the tt-diagrams for rrit = VJi-GeV/c^ 





with QCD (A) 


with QCD(B) 


Hww 


0.36 


0.35 


HwH 


0.77 


0.72 




0.34 


0.42 


sum 


1.47 


1.49 


V2 = 


0.57 


0.58 



Table 2: Corrections for the tt-diagrams for rrit — ITAGeV/c^. A) Results 
from full analytical formulae. B) Results after neglecting the evolution above 





without 


with QCD 
numerically 


with QCD 
negl. above mw 


Tja 

^WW 


1.91 


0.78 


0.78 


Ub 

^WW 


-0.041 


-0.006 


-0.007 


HwH 


0.39 


0.072 


0.074 


Hhh 


0.081 


0.012 


0.014 


sum 


2.34 


0.86 


0.86 


V3 = 




0.37 


0.37 



Table 3: Corrections for the ct-diagrams in units 10 ^ for = ITAGeV/c^ 



24 





For nil - 


= 150GV:l7c- 


For /;/( = 


= 200G'cl7r2 




without 


with QCD 


without 


with QCD 






numerically 




numerically 


Hww 


0.56 


0.34 


0.65 


0.40 


HwH 


1.00 


0.58 


1.50 


0.87 


Hww 


0.49 


0.24 


1.02 


0.51 


sum 


2.05 


1.16 


3.17 


1.78 


V2 = 




0.57 




0.56 



Table 4: The tt-case in units of 1 for values of rrit 





nit = 150GeV/c^ 


mt = 17 AGeV/c^ 


mt = 2mGeV/c^ 


nit = 50( 




without 


with QCD 


without 


with QCD 


without 


with QCD 


without 




Hww 

HwH 

Hhh 


0.56 
1.00 
0.49 


0.49 
0.84 
0.35 


0.61 
1.25 
0.72 


0.53 
1.04 
0.52 


0.65 
1.50 
1.02 


0.57 
1.26 
0.74 


0.87 
3.91 
8.55 




sum 

V2 


C 


1.82 


C 


1.81 


C 


.81 





Table 5: The tt-case for B-mesons in units of 1 
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